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Quantum ensembles form easily accessible architectures for studying various phenomena in quan¬ 
tum physics, quantum information science, and spectroscopy. Here we review some recent protocols 
for measurements in quantum ensembles by utilizing ancillary systems. We also illustrate these 
protocols experimentally via nuclear magnetic resonance techniques. In particular, we shall review 
noninvasive measurements, extracting expectation values of various operators, characterizations of 
quantum states, and quantum processes, and finally quantum noise engineering. 
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I. INTRODUCTION 

Unlike the classical measurements, measurements in 
quantum physics affect the dynamics of the system. 
Moreover, often a particular experimental technique may 
offer only a limited set of observables which can be di¬ 
rectly measured. The complete characterization of a 
quantum state or a quantum process requires, in general, 
a series of measurements of noncommuting observables - 
requiring repeated state preparation, and a large number 
independent measurements. In this article, we review re¬ 
cent progresses in the measurement of quantum ensem¬ 
bles and explain how to overcome the above challenges. 
Here we exploit the presence of an ancillary register in¬ 
teracting with the system that is to be measured. In the 
following section we show how to realize noninvasive mea¬ 
surements using ancillary qubits. Extracting expectation 
values of various types of operators and related applica¬ 
tions are described in section III. We then describe an 
efficient protocols for complete quantum state character¬ 
ization (section IV) and quantum process tomography 
(section V). We also narrate our experiments on noise 
engineering using ancillary qubits in section VI, and fi¬ 
nally we summarize all the topics in section VII. In all 
the sections, we illustrate the protocols experimentally 
using nuclear magnetic resonance (NMR) techniques. 


II. NONINVASIVE MEASUREMENTS 

A classical measurement can in principle be noninva¬ 
sive in the sense it has no effect on the dynamics of the 
system. The same is not true in general for a quantum 
system, wherein the process of measurement itself may 
affect the dynamics of the system. However, as explained 
below, ancillary qubits can be utilized to realize certain 
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FIG. I. Circuits describing invasive and non-invasive mea¬ 
surements. The first measurement is invasive in (a). The 
same is replaced by a noninvasive measurement in (b). The 
final measurements are invasive in both. 


quantum measurements without much disturbance, and 
hence extract probabilities or expectation values nonin- 
vasively to a great extent. Such measurements are often 
termed as noninvasive quantum measurements. 

Consider the example of a single qubit initialized in 
state p (Fig. 1 (a)) and a dichotomic observable Q with 
eigenvalues q G {0,1}. Suppose we need to extract the 
joint probabilities P{qi,q 2 ) at two time-instances, one 
before applying an unitary U and the other after apply¬ 
ing U (Fig. 1 (a)). To realize the first measurement 
noninvasively, we utilize an ancilla qubit (Fig. I (b)), a 
CNOT gate, and a final 2-qubit projective measurement in 
Q0Q basis. The CNOT operation copies the probabilities 
of P{qi) onto the ancilla qubit without projecting the 
system state. Denoting the first qubit as the ancilla and 
the second qubit as the system, the diagonal elements of 
the joint density operator pf store the joint probabilities, 
i.e., P(gi,g 2 ) = {QiQ 2 \pf\qiq 2 )- The joint probabilities 
can thus be extracted by a final strong measurement or 
by a diagonal density matrix tomography [1]. Knee et 
al [2] argued that since CNOT operator flips the ancilla 
qubit if the system qubit is in state |I), the circuit is 
not quite noninvasive. They also proposed a simple vari¬ 
ation, in which P(0, 0) and P(0,1) are measured using 
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CNOT operator, while P(l, 0) and P(l, 1) are measured us¬ 
ing an Anti-CNOT operator (which flips the ancilla only 
if the system qubit is in state |1)). In this procedure, 
called ideal negative-result measurement (INRM), all the 
joint probabilities are measured without flipping the an¬ 
cilla qubit, and therefore is considered more noninvasive. 
In the following we describe the application of INRM in 
studying Leggett-Garg inequalities. 

Leggett-Garg inequality (LGI) provides one way of 
distinguishing quantum behaviour from maerorealism 
Macrorealism is based on the following assumptions: (i) 
the object remains in one or the other of many possible 
states at all times, and (ii) the measurements are nonin¬ 
vasive, i.e., they reveal the state of the object without dis¬ 
turbing the object or its future dynamics. Leggett-Garg 
inequality (LGI) sets up macrorealistic bounds on linear 
combinations of two-time correlations of a dichotomic 
observable belonging to a single dynamical system [3]. 
Quantum systems are incompatible with these criteria 
and often violate bounds on correlations derived from 
them, thereby allowing us to distinguish the quantum 
behavior from maerorealism. Violations of LGI by quan¬ 
tum systems have been investigated and demonstrated 
experimentally in various systems [ 1 , 4-14]. An entropic 
formulation of LGI has also been introduced by Usha 
Devi et al. [15] in terms of classical Shannon entropies 
associated with classical correlations. We had reported 
an experimental demonstration of violation of entropic 
LGI (ELGI) in an ensemble of spin 1/2 nuclei using NMR 
techniques [ 1 ]. The simplest ELGI study involves three 
sets of two-time joint measurements of a dynamic ob¬ 
servable belonging to a ‘system’ qubit at time instants 
(^ 15 ^ 2)5 (^ 2 ,^ 3 ), and (^ 1 ,^ 3 ). The first measurement in 
each case must be noninvasive, and can be performed 
with the help of an ancilla qubit. 

Usha Devi et al. [15] have shown theoretically that 
for n-equidistant measurements on a spin-s system, the 
information deficit, [15] 
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FIG. 2. Information deficit D 3 versus d obtained using INRM 
procedure. The mean experimental D3 values are shown as 
symbols. The curves indicate theoretical D3. The horizontal 
lines at D3 = 0 indicate the lower bounds of the maerorealism 
territories (Figure reproduced from reference [1]). 


The resonance offset of ^^G was set to 100 Hz and that of 
^H to 0 Hz (on resonant). The two spins have an indirect 
spin-spin coupling constant J = 209.2 Hz. The NMR 
experiments were carried out at an ambient temperature 
of 300 K on a 500 MHz Bruker NMR spectrometer. 

As described in earlier, two sets of experiments were 
performed, one with CNOT and the other with anti-CNOT 
[1]. The joint entropies were calculated using the experi¬ 
mental probabilities and the information deficit (in bits) 
was calculated using the expression = 2 H{Q 2 \Qi) — 
H{Q 2 ,\Qi). The theoretical and experimental values of 
D 3 for various rotation angles 0 are shown in Fig. 2. 
According to quantum theory, a maximum violation of 
D 3 = —0.134 should occur at 6> = 7 r/ 4 . The correspond¬ 
ing experimental value, D 3 ( 7 r/ 4 ) = —0.114 ± 0.027, indi¬ 
cates a clear violation of ELGI. 

Our other experiments involving noninvasive measure¬ 
ments include (i) illustrating the inconsistency of quan¬ 
tum marginal probabilities with classical probability the¬ 
ory [ 1 ] and (ii) demonstrating that quantum joint prob¬ 
abilities can not be obtained from moment distribution 
[16]. 


III. EXTRACTING EXPECTATION VALUES 


where m G {I,-- - ,n}. The conditional probabilities 
in turn are calculated from the joint probabilities using 
Bayes theorem, 

P{(lj\(lk)P{qk) = P{qj,qk)‘ ( 2 ) 

We studied ELGI experimentally by treating the ^^G 
and ^H nuclear spins of ^^GHGl 3 (dissolved in GDGI 3 ) as 
the system and the ancilla qubits respectively (Eig. 2). 


Often experimental setups allow direct detection of 
only a limited set of observables and to extract their ex¬ 
pectation values. Eor example, in NMR only transverse 
magnetization operators {{cFx) and (cry)) are directly ob¬ 
servable via real and imaginary components of induced 
emf. Eig. 3 describes the circuits for measuring expecta¬ 
tion values of different types of operators using an ancilla 
qubit. Here the expectation values {ax) and (ay)) of an- 
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tronic states are 


FIG. 3. Circuits for (a) NMR measurements of {gx) and 
(cTj,), (b) standard Moussa protocol for expectation values of 
Hermitian-unitary operator A [17], and circuits for measur¬ 
ing expectation values of (c) a unitary operator U [18], (d) 
a projector P [18] (e) a diagonal Hermitian operator A, and 
(f) for measuring joint expectation values of noncommuting 
unitaries U and y[16]. 


cilia qubit reveal the expectation values of different types 
of operators acting on system qubit. Applications of such 
circuits are illustrated with the help of following exper¬ 
iments: (A) estimation of Franck-Condon factors, and 
(B) investigation of quantum contextuality. 


A. Estimation of Pranck-Condon factors 


^1 =pV2 + xV2 

'H2 = pV2 + (x - hf/2 + AE. 


(3) 


The FCF between \m) of the electronic ground state and 
|n) of the electronic excited state is given by [19], 




I i’mix)'lpn'{x,b)dx 


(4) 


where are the corresponding position 

wave-functions. 

Estimation of FCF, fm,n'^ is equivalent to measure¬ 
ment of expectation value of the projection = |m)(m| 
after preparing the system in excited state \n') since [18], 


fmx = {n'\m){m\n') 
{Pm)n'' 


(5) 


Three spin-1/2 ^^F nuclei of iodotrifluoroethylene 
(C 2 F 3 I) dissolved in acetone-De form a three-qubit NMR 
quantum simulator (see Fig. 4). Fi qubit is chosen to 
be the ancilla and the other two qubits chosen repre¬ 
senting the lowest four levels of the Harmonic oscillator 
[18]. The vibrational levels of the electronic ground state 
are encoded onto the spin states such that | 0 ) = jtt)? 
|1) = Iti), |2) = lit)^ and |3) = III). The preparation of 
excited state \n') can be achieved by first initializing the 
system in the corresponding state |n) of the electronic 
ground state and translating it in position from origin 
{x = 0) to the point x = b. This translation was be 
achieved by the unitary operator 


Urib) = e 




(6) 


Finally the expectation values (Pm) were measured ex¬ 
perimentally using the circuit shown in Fig. 3d, and then 
the FCFs fm,n were obtained using eqn. 5 [18]. The re¬ 
sults described in Fig. 4 display a good correspondence 
with the theoretically expected values indicating the suc¬ 
cess of the experimental protocols. 


Franck-Condon principle states that the transition 
probability between two vibronic levels depends on the 
overlap between the respective vibrational wavefunctions 
[19]. Franck-Condon factors (FCFs) dictate the intensi¬ 
ties of vibronic transitions, and therefore their estima¬ 
tion is an important task in understanding absorption 
and fluorescence spectra and related phenomena such as 
photo-induced dissociations [ 20 ]. 

We modelled the electronic ground and excited vi¬ 
brational levels as eigenstates of two harmonic poten¬ 
tials, Vi and V 2 respectively. To simulate the one¬ 
dimensional case, we choose the potentials Vi = x‘^/2 
and V 2 = {x — b)^/2 A- AE, which are identical up to an 
overall displacement b in position and/or in energy AE. 

Thus the vibrational Hamiltonians for the two elec- 


B. Investigation of quantum contextuality in a 
harmonic oscillator 

Quantum contextuality (QC) states that the outcome 
of the measurement depends not only on the system and 
the observable but also on the context of the measure¬ 
ment, i.e., on other compatible observables which are 
measured along with [ 21 ]. Consider the following NCHV 
inequality [ 22 ]: 

I = {AB eBC PCD- AD) 

= {AB) P {BC) P {CD) - {AD) < 2 (7) 

Hong-Yi Su et.al. [23] theoretically studied QC of 
eigenstates of ID-QHO. They introduced two sets of 
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FIG. 4. Experimental FCFs (circles) corresponding to 4- 
level harmonic oscillators versus the displacement b (in atomic 
units). The simulated FCFs (dotted lines) for the 4-level sys¬ 
tem and analytical FCFs (smooth lines) for infinite-level sys¬ 
tem are also shown for comparison. The dashed curve at 
the top of (b) corresponds to the normalization used. The 
thin vertical dashed lines at 6 = 2, 1 + mark the begin¬ 
ning of classically forbidden regions for /o,o', fo,i' respectively. 
Molecular structure and qubit labelling are shown in the inset 
(Figure reproduced from reference [18]). 


pseudo-spin operators, 

Tx = CFx0 1, Ty = (Jy, Tz = -(Jy C) (Jy, 

U = CTj; 0 az, Ty = 1 0 Gy, U = ® (8) 

where, 1 is 2 x 2 Identity matrix. Using these operators 
they defined the observables, 

^ ~ B = ~ C) {Ci^dz S^CTx) 

C = r^;, D = CyY T SyY^ = (7x 0 {CyCTz (9) 

where = cos/3, Sj^ = sin/3, Cy = cos 77 , and Sy = sin?]. 
Here operators M, H, (7, D are unitary & Hermitian and 
accordingly have eigenvalues ± 1 , with (B^C)^ 

{C,D) and {D^A) forming compatible pairs. Hong-Yi 
Su et.al. have shown that. 




(-7r/4, -37r/4)o 
(37r/4,7r/4)i 
(7r/4,37r/4)2 ^ ’ 

(-37r/4, -7r/4)3 


where, is the expression on LHS of inequality 

7, I = 0,1,2 and 3, and, \0)qho,\^)qho,\‘^)qho and 
\^)qho are first four energy eigenstates of ID-QHO. 

We encoded the first four energy eigenstates of ID- 
QHO onto the four Zeeman eigenstates of a pair of spin- 
1/2 nuclei. The circuit shown in Fig. 3f was used to ex¬ 
tract the expectation value of observables (HH, HC, CD^ 
and DA) in a joint measurement. We used three ^^F nu¬ 



FIG. 5. lo, Ii, I2, and I3 represent evaluation of inequality 
7 for eigenstates \0)qho, | l ) QifO , 1^2)qho, and \3)qho respec¬ 
tively [24]. The curved surface represents theoretical values, 
and the points are experimental values. Flat planes at / = 2 
and / = — 2 represent classical bounds. 


clear spins of trifluoroiodoethylene dissolved in acetone- 
D 6 (see inset of Fig. 4) as the 3-qubit register. The first 
spin, Fi, was used as an ancilla qubit, and other spins, F 2 
and F 3 , as the system qubits. The results are shown in 
Fig. 5. The maximum theoretical violation is 2\^ = 2.82 
[24]. The experimental value of maximum violation for 
Io,Ii,l2, and I3 are 2.40 ± 0.017, 2.45 ± 0.025, 2.39 ± 
0.016, and 2.42 ±0.026 respectively [24]. There is a clear 
violation of the classical bound. Reduced violation than 
the theoretical value is due to T 2 decay and inhomogene¬ 
ity in Radio Frequency (RF) pulses. 


IV. ANCILLA ASSISTED QUANTUM STATE 
TOMOGRAPHY 

In experimental quantum information studies. Quan¬ 
tum State Tomography (QST) is an important tool that 
is routinely used to characterize an instantaneous quan¬ 
tum state [25]. QST can be performed by a series of mea¬ 
surements of noncommuting observables which together 
enables one to reconstruct the complete complex den¬ 
sity matrix [25]. In the standard method, the required 
number of independent experiments grows exponentially 
with the number of input qubits [26, 27]. Anil Kumar 
and co-workers have illustrated QST using a single two- 
dimensional NMR spectrum [28]. Later Nieuwenhuizen 
and co-workers showed how to reduce the number of in¬ 
dependent experiments in the presence of an ancilla reg¬ 
ister [29] . We referred to this method as Ancilla Assisted 
QST (AAQST) and experimentally demonstrated it us¬ 
ing NMR systems [30, 31]. AAQST also allows single 
shot mapping of density matrix which not only reduces 
the experimental time, but also alleviates the need to 
prepare the target state repeatedly [31]. 

To see how AAQST works, consider an input register 
of ?r-qubits associated with an ancilla register consisting 
of ft qubits. The dimension of the combined system of 
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FIG. 6. Minimum number of independent experiments re¬ 
quired for QST (with zero ancilla) and AAQST (Figure re¬ 
produced from reference [31]). 

n = n + n qubits is TV = NN^ where N = 2^. A com¬ 
pletely resolved NMR spectrum yields hN real param¬ 
eters. We assume that the ancilla register begins with 
the maximally mixed initial state, with no contribution 
to the spectral lines from it. The deviation density ma¬ 
trix of the combined system is p = p ( 8 ) 1/TV. To perform 
AAQST, we apply a non-local unitary of the form, 

N-l 

Uk = VY,Uka®\a){a\. ( 11 ) 

a=0 

Here Uka is the kth unitary on the input register depen¬ 
dent on the ancilla state \a) and V is the local unitary 
on the ancilla. The combined state evolves to 

= UkpUl 

= jY. Pmm-Uka\m){m'\Ul^®V\a){a\V^. (12) 

m,m' ,a 

Intensity of NMR spectrum is proportional to the observ- 

n 

able ^jx + The spectrum of the combined sys- 

T=i 

tern yields fiN linear equations. The minimum number 
of independent experiments needed is now 0{N‘^/{fiN)). 
Choosing N ^ N, AAQST needs fewer than 0(TV/n) ex¬ 
periments required in the standard QST. In particular, 
when fiN > TV^, a single optimized unitary suffices for 
QST. Fig. 6 illustrates the minimum number {K) of ex¬ 
periments required for various sizes of input and ancilla 
registers. As illustrated, QST can be achieved with only 
one experiment, if an ancilla of sufficient size is provided 
along with [31]. 


To demonstrate this procedure experimentally, we used 
three ^^F nuclei and two nuclei of l-bromo-2,4,5- 
trifluorobenzene (BTFBz) partially oriented in a liquid 
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FIG. 7. AAQST results for the state described in the text [31]. 
The reference spectrum is in the top trace. The real (middle 
trace) and the imaginary spectra (bottom trace) are obtained 
in a single shot AAQST experiment. The bar plots corre¬ 
spond to theoretically expected density matrices (top row) 
and those obtained from AAQST experiments (bottom row). 
The molecular structure and the labelling scheme of BTFBz is 
shown in the center (Figure reproduced from reference [31]). 


crystal namely, N-(4-methoxybenzaldehyde)-4- butylani- 
line (MBBA) (Fig. 7) [31]. We chose the three ^^F nuclei 
forming the input register and two nuclei forming the 
ancilla register. 

Fig 7 shows the experimental results corresponding to 
a particular density matrix obtained by applying unitary 
^0 = To {^Yy , with To = 2.5 ms, on ther- 

mal equilibrium state. The real and imaginary parts of 
the reconstructed density matrix along with the theo¬ 
retically expected matrices are shown below the spectra 
in Fig. 7. Fidelity of the experimental state with the 
theoretical state was 0.95. The entire three-qubit den¬ 
sity matrix with 63 unknowns was estimated by a single 
NMR experiment [31]. 


V. SINGLE SCAN PROCESS TOMOGRAPHY 

Often one needs to characterize the overall process 
acting on a quantum system. Such a characterization, 
achieved by a procedure called quantum process tomog- 
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raphy (QPT), is crucial in designing fault-tolerant quan¬ 
tum processors [32, 33]. QPT is realized by considering 
the quantum process as a map from a complete set of 
initial states to final states, and experimentally charac¬ 
terizing each of the final states using QST [26]. Since 
QST by itself involves repeated preparations of a target 
state, QPT in general requires a number of independent 
experiments. Therefore the total number of independent 
measurements required for QPT increases exponentially 
with the size of the system undergoing the process. 

The physical realization of QPT has been demon¬ 
strated on various experimental setups [34-48]. Several 
developments in the methodology of QPT have also been 
reported [49, 50]. In particular, it has been shown that 
ancilla assisted process tomography (AAPT) can charac¬ 
terize a process with a single QST [36, 37, 51, 52]. By 
combining AAQST and AAPT, we showed that entire 
QPT can be carried out with a single ensemble measure¬ 
ment [53]. We referred to this procedure as ‘single-scan 
quantum process tomography’ (SSPT) [53]. 

In the normal QPT procedure, the outcome of the pro¬ 
cess 5 is expanded in a complete basis of linearly indepen¬ 
dent elements {pi, P 25 * * * 5 Pn^} well as using operator- 
sum representation, i.e., 

^{Pj) ~ ^ V ^jkPk — ^ ^ ^iP^i ' (13) 

k i 

The complex coefficients Xjk can be extracted using QST. 
We can utilize a fixed set of basis operators {Em}, and 
express Ei = CimEm so that 

^{p) — ^ EmpEj^Xmn, 

mn 

where Xmn = ^im^in' The y matrix completely char¬ 
acterizes the process 5 . Since the set {pk} forms a com¬ 
plete basis, it is also possible to express 

k 

where can be calculated theoretically. Using eqns. 
13-15 and using the linear independence of {pk}, we ob¬ 
tain 


Px = A, (16) 

from which y-matrix can be extracted by standard meth¬ 
ods in linear algebra. 

A comparison of QPT, AAPT, and SSPT procedures 
for a single qubit process is presented in Fig. 8 [53]. 
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FIG. 8. Illustrating (a) single-qubit QPT requiring a total of 8 
NMR scans, (b) AAPT requiring 2 NMR scans, and (c) SSPT 
requiring a single-scan NMR experiment [53]. In each case, 
dotted lines are used to indicate the single-quantum elements 
of the density matrix which are directly observable. Other el¬ 
ements are observed by converting them to observable single¬ 
quantum coherences by using certain unitary operations in a 
subsequent scan(s) (Figure reproduced from reference [53]). 


n 

TTqpt 

ATaapt 

(«a) 

Msspt 

{nA,nB) 

1 

8 

2 

(1) 

1 

(1,1) 

2 

32 

4 

(2) 

1 

(2,2) 

3 

192 

11 

(3) 

1 

(3,3) 

4 

1024 

32 

(4) 

1 

(4,5) 

5 

7168 

103 

(5) 

1 

(5,6) 


TABLE I. Comparison of number of scans and number of an¬ 
cilla qubits (in parenthesis) required for n-qubit QPT, AAPT, 
and SSPT. 


experimentally obtained y matrices for certain quantum 
processes using the single scan procedure are shown in 
Fig. 9. 


Estimates of number of measurements for a small number 
of qubits shown in the first column of Table 1 illustrate 
the exponential increase of Mqpt with n. 

The experimental demonstration of a single-qubit 
SSPT was carried out using iodotrifiuoroethylene dis¬ 
solved in acetone-De as a 3-qubit system Fig. 9 [53]. The 


VI. ANCILLA ASSISTED NOISE 
ENGINEERING 

Preserving coherence is a very important aspect to re¬ 
alize quantum processors, and hence various techniques 
have been developed to suppress decoherence. They 
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FIG. 9. The barplots showing experimental X"iiicLtrices for 
various quantum processes obtained using SSPT. In each case, 
the left and right barplots correspond to the real and imag¬ 
inary parts respectively, and the hdelities are indicated in 
parenthesis (Figure reproduced from reference [53]). 


include dynamical decoupling (DD) techniques [54-57], 
quantum error correction [58], adiabatic quantum com¬ 
putation [59] , and use of decoherence-free subspaces [60] . 
Earlier Teklemariam et al. introduced artificial decoher¬ 
ence by achieving irreversible phase damping via con¬ 
stant perturbation of the environment qubits (by random 
classical fields), and thus mimicking a large dimensional 
environmental bath. Such experiments provide insights 
about decoherence processes and may pave the way for 
improving decoherence suppression techniques. 

In our work we simulated such a decoherence process 
on a NMR simulator with two qubits, where one acts as 
the system and the other as environment [61]. We then 
subjected the system qubit to certain DD sequences and 
observed their competition with the engineered decoher¬ 
ence through noise spectroscopy. 

The two qubit register was initially in the product 
state, p(0) = p^(0)(8)p®(0). Here p^(0) is the system state 
and p^(0) is the environment state. We chose and 
nuclear spins in ^^C-labelled chloroform (^^CHCls dis¬ 
solved in CDCI 3 ) as the system and environment qubit 
respectively. The NMR Hamiltonian is 

n = (17) 

where Ug and Ue are the resonant frequencies of the sys- 



FIG. 10. Logarithm of transverse magnetization log(Mcc) as a 
function of time under different pulse sequences as indicated 
in the legend. Here r = 3.2 ms, F = 25 kicks/ms, and e G 
[0°, 1°]. The T2 values for various pulse sequences are shown 
in the legend (Figure reproduced from reference [61]). 


tern and the environment qubits respectively, J is the 
coupling strength between the two, and cr|, are the 
Pauli operators. In a total duration T, the propagator 
U = entangles the system qubit with the environ¬ 

ment qubit via the interaction J. We engineered decoher¬ 
ence by a series of RF kicks of arbitrary angles e G [—0^0] 
on the environment qubit. These kicks induced artifi¬ 
cial decoherence on the system qubit. Teklemariam et 
al. proved that induced decoherence of the system qubit 
depends on the kick-rate F, range of kick-angles and 
coupling strength J [62]. Their model predicted that for 
small kick-angles e and for lower kick rates F, decoher¬ 
ence rate I/T 2 increases linearly with F. After a certain 
value of F, I/T 2 saturates, and then onwards, it decreases 
exponentially with F. 

Our experimental results for e G [0°,1°] and F = 25 
kicks/ms are shown in Fig. 10 (indicated by stars). For 
comparison we have also shown the decay of magnetiza¬ 
tion without kicks (indicated by filled circles). Evidently, 
the kicks on environment have introduced additional de¬ 
coherence thereby increasing the decay of system coher¬ 
ence [61]. 

Dynamical decoupling attempts to inhibit decoherence 
of system by rapid modulation of the system state so as 
to cancel the system-environment joint evolutions. The 
two standard DD sequences are: (i) CPMG [54] and (ii) 
UDD [57]. CPMG consists of a series of equidistant tt 
pulses applied in the system qubits. In an N-pulse UDD 
of cycle time tc, the time instant tj of 7r-pulse is 


given by tj = tc sin' 
ments for t, 




. The results of the experi- 


2(Ar+i)_ 

22.4 ms and with different kick-parameters 
are shown in Eig. 10. The competition between kicks- 
induced decoherence and DD sequences can be readily 
observed [61]. 

Noise spectroscopy provides information about noise 
spectral density, which is the frequency distribution of 






















VII. SUMMARY 
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FIG. 11. Experimental spectral density profiles with different 
kick-angles (as indicated in the legend) and with kick-rates 50 
kicks/ms (top trace) and 25 kicks/ms (bottom trace). In both 
the traces, experimental spectral density profile without kicks 
is also shown for comparison. The smooth lines correspond 
to fits with one or two Gaussians (Figure reproduced from 
reference [61]). 

noise and is helpful in not only understanding the perfor¬ 
mance of standard DD sequences, but also in optimizing 
them [63-65]. In the limit of a large number of tt pulses, 
the CPMG filter function resembles a delta peak at cj, 
and samples this particular spectral frequency. The am¬ 
plitude of the noise S{uj) can be determined by using the 
relation [66] S{uj) 47f(T)- Thus by measuring T2{uj) 
for a range of cj = tt/t values, we can scan the profile of 
S{w). 

The experimental spectral density profiles of only natu¬ 
ral decoherence (lowest curve in each sub-plot), and with 
kicks of different kick-parameters are shown in Fig. 11 
[61]. Clearly the effect of kicks is to increase the area 
under the spectral density profiles and thereby leading 
to faster decoherence. Moreover, for a given kick-rate 
F, larger the range of kick-angles, higher is the spectral 
density profile [61]. 


In this review article, we described several recent pro¬ 
tocols for efficient measurements on quantum systems, 
and illustrated their NMR implementations. 

In section II, we described ancilla-assisted noninvasive 
measurements, where the measurement result of an inter¬ 
mediate observable was temporarily stored in an ancilla 
qubit. A final joint-measurement of the system-ancilla 
register revealed the joint probabilities in a noninvasive 
way. We also showed the application of this technique in 
studying entropic Leggett-Garg inequality [1]. 

In section III, we described extracting expectation 
values of various types of operators. Applications of 
these methods are illustrated in the estimation of Franck- 
Condon coefficients and in the investigation of quantum 
contextuality [18, 24]. 

In section IV and V we described efficient ways to char¬ 
acterize quantum states and quantum process by exploit¬ 
ing ancilla qubits. We also illustrated single-scan quan¬ 
tum state tomography as well as single-scan quantum 
process tomography using NMR systems. These tech¬ 
niques not only alleviate the need of repeated measure¬ 
ments, but also allow the study of random states or dy¬ 
namic processes [31, 53]. 

Finally in section VI we described ancilla-assisted noise 
engineering, where random fields applied to the an¬ 
cilla qubits cause controllable decoherence on the system 
qubits. We illustrated this phenomena using a two-qubit 
NMR system, and studied the engineered decoherence by 
measuring noise spectrum [61]. 

Although we have used NMR techniques to demon¬ 
strate the above protocols experimentally, these proce¬ 
dures are quite general in nature, and can easily be 
adopted to other experimental techniques as well. 
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